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Abstract. The electrical and optical properties of ordered passive arrays, constituted of inductive and 
capacitive components, are usually deduced from Kirchhoff's rules. Under the assumption of periodic 
boundary conditions, comparable results may be obtained via an approach employing transfer matrices. 
In particular, resonances in the dielectric spectrum are demonstrated to occur if all eigenvalues of the 
transfer matrix of the entire array are unity. The latter condition, which is shown to be equivalent to the 
habitual definition of a resonance in impedance for an array between electrodes, allows for a convenient 
and accurate determination of the resonance frequencies, and may thus be used as a tool for the design 
of materials with a specific dielectric response. For the opposite case of linear arrays in a large network, 
where periodic boundary condition do not apply, several asymptotic properties are derived. Throughout 
the article, the derived analytic results are compared to numerical models, based on either Exact Numerical 
Renormalisation or the spectral method. 



PACS. 77.22.-d Dielectric properties of solids and liquids - 78.20.-e Optical properties of bulk mate- 
rials and thin films - 78.20.Bh Theory, models, and numerical simulation - 41.20.-q Applied classical 
electromagnetism 



' 1 Introduction 

' Recent experimental advances have put to living the 
, speculation about materials with negative refraction in- 
' dex, initiated in 1968 by Victor Veselago on purely theo- 

■ retical grounds 2^. Such materials, starting to be at reach 
nowadays primarily for the microwave region, manifest ex- 

' citing and unconventional phenomena ranging from the 
, Inverse Doppler effect jS, to the possibility of diffraction- 

■ free imaging 

' Were such "left-handed" materials, as Veselago called 
I them, available for any frequency domain, they would un- 
' doubtedly revolutionise optics. However, a negative refrac- 
, tion index requires simultaneously a negative permittivity 
' and a negative permeability, implying that the material 
has dielectric and magnetic resonances in the same fre- 
quency domain — a property which turns out to be ex- 
tremely rare and which is partly responsible why such 
materials have remained undiscovered for almost three 
decades. 

Unlike to the optical domain, where the electric and 
magnetic resonances are generally separated in frequency 
by several orders of magnitude, negatively refracting ma- 
terials have recently been engineered for narrow bands in 
the microwave domain. One of the main ingredients of 
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these "metamaterials" are regular arrays of so-called split- 
ring resonators 12] . 

In this paper, we will focus on the dielectric part of 
the response of arrays of similar resonators in square lat- 
tices. One goal is to provide tools for the location of the 
resonances in the spectrum, a task which cannot simply 
be performed by looking at the building blocks of such 
an array IP. At the same time, the theoretical approaches 
presented in this paper allow for a deeper understanding 
of the resonance spectrum — certainly an advantage in 
the quest of materials with a tailored electromagnetic re- 
sponse. Throughout the paper, the analytical results are 
compared to numerically calculated spectra, obtained ei- 
ther via Exact Numerical Renormalisation (ENR) 
|H] or from the spectral method |51llUllll| formulated with 
Green's functions pitTMl] . 

An archetype of the circuits to be dealt with in this 
paper is depicted in Fig. This "cluster" consists of two 
coils in series (each of inductance L and with no ohmic 
resistance, R — Q) embedded in the simplest imaginable 
network made of three capacitors (each of capacity C). 
As will be seen in the following, this circuit contains — 
despite its simplicity — many of the essential features of 
the more complicated arrays to be studied in the following 
sections. 

With no external voltage applied to the plates, the 
state of the system can be obtained by straightforward 
application of Kirchhoff's rules: the loop rule yields a dif- 
ferential equation for the current through any of the coils 
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Fig. 1. One of the simplest non-trivial clusters, consisting of 
two self-induction coils embedded in a capacitor network. 



Ik in terms of the currents flowing through the vertical and 
the horizontal capacitor belonging to the same loop; the 
latter currents are subsequently eliminated by KirchhofF's 
junction rule, resulting in a homogeneous system of two 
ordinary differential equations for the currents through the 
coils. Applying a time-dependent external voltage Vc->it(t) 
to the plates introduces an inhomogeneity in the equa- 
tions, which finally read 



dt2 



uj; 



-2 1 
1 -2 



(1) 



with u)o — 1/ \J LC and /oxt (i) the total net current fiowing 
through the sample. 

Eq. can be solved by standard means yielding the 
homogeneous eigenmodes 



ihW 



Cl 



C2 



(2) 



where c\ and C2 are arbitrary complex amplitudes. The 
application of a sinusoidal external voltage of frequency 
i7, implying an external current /cxt(i) = -^oxt exp(if2t), 
allows for one additional solution 



linh(^) — -^c: 




for i7 7^ tjQ, 
for fl = u}Q. 



(3) 



Four observations can be made at this point: (i) in 
networks consisting of two species of components only, 
the state of the system is fully described by a set of dif- 
ferential equations for the currents flowing through the 
minority components (in our case the coils). The con- 
tribution of the majority components (here, the capaci- 
tors) to the Kirchhoff rules is purely algebraic, and may 
be eliminated from the system of equations, (ii) The ex- 
ternal voltage can only excite the first resonance, of fre- 
quency • The second resonance, -v/Si-^O: is antisymmet- 
ric {1\ = —I2) and thus orthogonal to the currents which 
are always induced symmetrically in the present config- 
uration of the voltage plates, (iii) The applied voltage 
equals the sum of the two voltage drops in the coils, i.e. 
V^cxt = L'Sti^^ + ^2)- Off resonance, this yields a total 
impedance of Z = K:xt/^oxt = 2iLnu}Q/ {uJq — H'^) between 
the plates. This result may be generalised: as a resonance 
is approached, the internal currents of the sample diverge 
for finite applied voltage. This may be associated with 
an infinite conductivity of the sample's (internal) compo- 
nents. At the same time, however, a divergent impedance 



and thus zero conductivity is measured between the plates, 
(iv) In reality, the coils have a small but finite ohmic resis- 
tance R. The eigenfrequencies are thus shifted slightly in 
the complex plane towards small positive imaginary parts 
which damp out the sample's free modes, eq. In ad- 
dition to its imaginary part, the impedance Z acquires a 
real part consisting of a narrow lorentzian peak centred at 
the resonance frequency. 

The paper is organised as follows: in Sec. [2 the res- 
onance spectra of regular arrays with not too compli- 
cated unit cells are calculated directly from the solution 
of Kirchhoff's rules. In Sec. |31 more complicated arrays 
are tackled with an approach based on transfer matrices. 
Regular one-dimensional (ID) arrays in an infinite two- 
dimensional (2D) lattice are examined in Sec.^ and the 
results compared to the formerly discussed 2D clusters. 
Finally, the Appendix ^ is devoted to a simple physical 
approximation, based on a dipole scenario in a 2D environ- 
ment, which turns out to be helpful for the interpretation 
of the spectra of linear clusters. 



2 Direct solution of Kirchhoff's rules 

In this section, we will tackle simple regular ID and 2D 
binary arrays by a direct solution of Kirchhoff's rules. The 
results are then to be compared to those obtained by Ex- 
act Numerical Renormalisation (ENR) 6,7,8 . The strat- 
egy of the latter algorithm consists of eliminating the net- 
work sites one by one, while renormalising the impedance 
between all couples of former neighbours of the eliminated 
site such that the global impedance remains invariant, un- 
til the electrodes are connected by just one bond, which 
then carries the whole network's impedance. 




Fig. 2. A ladder-shaped circuit containing A'^ coils in series 
on one leg, and the same number of bonds with no impedance 
on the other. AH N + 1 steps of the ladder are capacitors. 



To begin with, we are going to summarise the solution 
of Kirchhoff's rules for the ID array shown in Fig. 12 This 
ladder-shaped circuit consists of a horizontal line of N 
coils, each of inductance L, which is connected by -I- 1 
capacitors, each of capacity C, to another horizontal line 
with no resistance at all. The loop rule states that the 
voltage drop around a closed loop is zero; applying it to 
the fc.th mesh and deriving with respect to time yields 
d^/fc/dt^ = ujQ(ii. — ii._i), withwo = I/VLC. The currents 
through the capacitors, ik, can be eliminated using the 
junction rule ik + Ik = h+i, where Ik is assumed zero if 
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k is out of range. In total, we get 



(4) 



where I — (Ii, . . . , In) is the coil current vector, and D 
the ID lattice Laplacian, i.e. a tridiagonal matrix with 
entries —2 for all diagonal elements and 1 for all non-zero 
off-diagonal elements. 

The usual ansatz I exp [icut) converts the differential 
equation Q into an eigenvalue problem for the tridiagonal 
matrix D-t-cj^l, where uj = uj/ujq stands for the frequency 
in units of luq. Its determinant can be calculated explic- 
itly: for a A'^-dimensional tridiagonal matrix M.]^{x) with 
diagonal elements 2x, and first sub- and super-diagonal el- 
ements 1, expansion of the determinant by minors yields 



Dn{x) = 2xDN-iix) - Dn^2{x) 



(5) 



where Dn{x) = detMjv(a;). We recognise in ||SJ) the recur- 
rence relation generating the Chebyshev polynomials. The 
required initial conditions, Dq{x) = 1 and Di(x) = 2x, 
ties us down to the Chebyshev polynomials of the second 
kind: 

, / .^ sin[(iV + 1)6'1 

Un{x ^co&e) - . ' (6) 



The roots of Un {cos 9), which occur at 9„ 



""^ deter- 

mine the eigenfrequencies. In the present example, where 
2x = 2 cos 9m = —2 + oj'^, after renumbering the eigen- 
states according to their frequency, we have 



= ±2 sin 



2(iV+ 1) 



(for m = l...N). 



(7) 



For infinite ladder length, the dispersion relation uj{k) ver- 
sus wave number k — is shown in the first graph of 
Fig. The associated density of states 



lim 

N- 







(8) 



shown in the second graph, is finite for — 0, since 
the sample's lowest resonance frequencies go to zero for 
N ^ oo. The fact that p{uj) displays Van Hove singu- 
larities for (D = ±2 should not be misinterpreted in the 
sense that the sample mainly resonates in these frequen- 
cies when exposed to a multi- frequency input: on the con- 
trary, as will be pointed out in Sec. O for reasons of 
symmetry there is very little overlap between the volt- 
age applied to the electrodes, generating an overall cur- 
rent with symmetry loxt ~ (1; I7 • ■ • 1 1)"^: find the highest 
excited states, with m close to N. 

In the following, we are going to study regular 2D ar- 
rays generated from a unit cell which tiles the surface be- 
tween the electrodes. At first, the unit cell is replicated N 
times in x-direction, i.e. from the left towards the right 
electrode. The resulting horizontal ladder is then stacked 
in 2/-direction until the width of the electrodes is fully 
covered, assuming periodic boundary conditions. In other 
words, the 2D arrays to be studied are cylinders of length 




1 






1 



Fig. 4. Staircase-array obtained by replicating the unit cell of 
Fig- Eld) iV = 3 times in a;-direction. The resulting structure is 
replicated in j/-direction until the electrodes are fully covered. 
Periodic boundary conditions are assumed in y-direction. 



(a) 



1 



(b) 



(0 j_ j_ 




(A) 





Fig. 5. Unit cells for several 2D arrays (left column) and 
corresponding ID building blocks (right column). 



N (times the length of the unit cell). At either end, they 
are connected to an electrode ring via at least one rank 
of pure capacitor bonds, which avoids short circuiting and 
keeps the influence of the electrodes on the dielectric spec- 
tra as low as possible. 

The building block which generates an array resem- 
bling the ladder in Fig. |21 is shown in Fig. |Sfa). Similar 
arrays can be generated from the other unit cells in Fig.|21 
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Fig. 3. Dispersion relations Lu{k) (left column) and corresponding densities of states p{uj) (right column). The two graphs in 
the first line correspond to the ladder of Fig. |2| (continuous lines — eqs. 10 and lO) and to its 2D counterpart shown in Fig.|3a) 
(dashed lines). The two graphs in the second line where obtained for the arrays of Fig. |^b) (continuous lines — eqs. Q and 
TM). and Fig.l^c) (dotted lines — eq. lO). 



the array shown in Fig. ^ was obtained from the unit cell 
represented in Fig.E^d) with = 3. 

As illustrated in the right column of Fig. these 2D 
arrays can be mapped onto ID ladders: translational in- 
variance along the electrodes in steps of an integer times 
the height of the unit cell allows to "bend down" the un- 
bound components in the top row and to attach them to 
the corresponding site at the bottom. Therefore, the only 
difference between the original ladder of Fig. |21 and the 
ladder in Fig. EJa) resumes to doubling the vertical capac- 
itors, and thus dividing the frequency range by a factor 
of y/2 in eqs. iQ and ©. (The corresponding dispersion 
relation and density of states are plotted with dashed lines 
in the first row of Fig. |3Jl . 

The array in Fig.jSfb) can be tackled analogously. The 
additional capacitors only contribute to the diagonal ele- 
ments of the tridiagonal matrix, and instead of eq. I@J one 

2 

obtains d^I/di^ = ^[D — 2 • The resulting resonance 
frequencies. 



<i„ = ±^l + 2sin2^^^^-^ (formal... TV), (9) 

are plotted, for = 5, with dotted lines in Fig. El A 
corresponding ENR calculation for the total conductiv- 
ity of the array between electrodes (continuous lines in 



Fig. El) detects the resonance frequencies with m even, as- 
sociated with symmetric eigenmodes. The antisymmetric 
eigenmodes are orthogonal to the current vector induced 
by the plates, and remain hence invisible in this calcula- 
tion. 

The dispersion relation obtained from eq. lO for N 
oo, albeit with k — finite, along with the correspond- 
ing density of states. 



is illustrated by the continuous lines in the graphs in the 
second row of Fig. 13 This time, even for infinite chain 
length, there are no low-lying resonances, and a gap spreads 
between —1 and -1-1. 

The last array to be treated with this method can be 
obtained from a pure capacitor lattice by replacing ev- 
ery second horizontal capacitor on every second line by a 
coil. In this array, which corresponds to the unit cell in 
Fig. ISIc), meshes with one coil and three capacitors al- 
ternate horizontally with pure capacitor meshes. For the 
latter, Kirchhoff's rules reduce to purely algebraic rela- 
tions; these can be removed from the system of differential 
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Fig. 6. Impedance between the electrodes for an array gen- 
erated from the unit cell in Fig. Elb) . The sample consists of 6 
lines, of A'^ = 5 inductors in series each, substituting the central 
piece of every second row in a 12 by 12 capacitor lattice. The 
real part of the ENR impedance, calculated for L = C = 1 and 
Rcoa = 10"^, is plotted with continuous lines. The dotted lines 
indicate the resonance frequencies obtained analytically from 
eq. ||^ for Ai' = 5 and m = 1 . . . 5. 



equations, and one ends up with 
/-23 1 
1 -22 1 



dt^ 



"^0 

12 



\ 





V 



1 -22 1 
1-22/ 



(11) 



where I is the current vector through the N coils of the 
cluster. The matrix in eq. (|ll|l differs slightly from the ha- 
bitual form, since not all diagonal elements are the same. 
The asymmetry in the upper left corner translates the fact 
that the corresponding ladder (right column of Fig. Oc)) 
starts with a loop of 3 capacitors and 1 coil, but ends with 
a 4-capacitor loop. 

The solution of eq. (|ll|l is analogous to the generic 
case, eq. (@}, albeit leading to a slightly more complicated 
recursion relation. A cumbersome but straightforward cal- 
culation gives the dispersion relation 



/ 5 1.9 riv 

±\ - + - sm^ 

^'3 3 2N- 



(form = l...iV), (12) 



with an associated density of states showing narrow bands 
in the range of PI = ^5/3 ~ 1.29 to V2 ~ 1.41. As can 
be seen from the dotted curves in the last two graphs of 
Fig. 13 the dispersion flattens substantially, and the bands 
in the density of states are squeezed with respect to the 
former case, without pure capacitor meshes (continuous 
lines in the same graphs, cf. eqs. (O and (|10|l 'l. 

This behaviour translates the tendency of the reso- 
nances to localise on individual loops as more and more 
pure capacitor loops are inserted in the ladder. In the limit 



of infinite distance between coils, 
coupled LC circuits of the type 



one is left with N de- 



C. 



inf 



c 



(13) 



with 



c c 



c c c = c c 

V — I — V — I — V — I — V — I — — I — in 

c c c 



(14) 



The solution of ^ is dnf = ^Cv{l + + 2C/Cv) 
which in our case, where Cy = 2C, reduces to Cinf = 
C(l + \/2). All N LC circuits of type lO resonate thus 

at the same Hi — ±^2^2 — 1 ~ ±1.35, and the density 
of states of the system reduces to delta functions at these 
frequencies which lie in the centre of the narrow dotted 
bands shown in the fourth graph of Fig. |3| 

The main advantage of the direct solution of Kirch- 
hoff's rules, presented for several arrays in this section, 
are its analytical results and the physical insight it pro- 
vides into the structure of the resonance spectra. On the 
other hand, for increasingly complex arrays, the method 
requires — if feasible at all — more and more cumbersome 
calculations for the solution of the recursion relations it 
relies on. We will therefore present an algorithm which 
circumvents this problem in the next section. 



3 Transfer matrix method 

The analysis of the resonance spectra of the arrays dis- 
cussed so far — which due to their translational invari- 
ance along the electrodes reduce to effective ID problems 
— can be rephrased very efficiently using transfer matri- 
ces. 

To illustrate the concept of a transfer matrix, consider 
the quadrupole 



1 



(15) 

Its incoming and outgoing currents and voltages are con- 
nected to each other via 



TJi 



(16) 



where = {Uk,Ik, Uk', Ik') is the vector of state at each 
end k = 1,2, i.e. a combination of currents / and voltages 
U = Vao (the latter, for convenience, multiplied by the 
conductivity ctq of the lattice's majority component); T 
represents the dimensionless transfer matrix. 
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For quadrupoles exclusively made of passive elements, 
there are two conservation laws. The first states that the 
potentials are determined up to an arbitrary additive con- 
stant Uq. Increasing all incoming potentials by Uq shifts 
the outgoing ones by the same amount: hence (1, 0, 1, 0)"^ 
is a right eigenvector of the generally non-symmetric ma- 
trix T, with corresponding eigenvalue /i = 1. Secondly, 
the continuity equation states that the total incoming and 
outgoing current have to be the same: hence, (0, 1, 0, 1) is 
a left eigenvector of T, again associated with /i = 1. 

These conservation laws may be easily verified since all 
passive quadrupoles can be assembled by combining two 
prototypes. 



(A) 



a 



(17a) 



(B) 



with corresponding transfer matrices 



10 

1-i 



Vo 1 / 















and Tb = 




1 




V 

1 


i) 


/ 


yl 





-?7 





(17b) 



where rj — a/ao (and 77' = a'/ao) denote the conductivi- 
ties of the components (in units of (Tq). Furthermore, since 
the transfer matrix T of any passive quadrupole can be 
written as a product of matrices of type Ta and Tb, one 
always has det T = 1 — a relation that will be of some 
importance in the following. 

Naively one could define a resonance as "if what goes 
in, comes out" . This amounts to resolving 



AU 



(18) 



where (arbitrarily) leg 2 has been grounded, U2 = 0; AU 
denotes the voltage drop between the two sides of the 
quadrupole, and w stands for the potential difference be- 
tween two legs on the same side of the quadrupole. When 
placed between electrodes, a resonance would be seen as 
an impedance pole, implying that at finite AU the net 
current through the sample is zero, Ii = I2 + I2' = 0. 
Apart from the technical difficulty that H18|l relies on the 
solution of a singular system of equations, this method 
nicely yields the lowest eigenfrequency of the quadrupole. 

By contrast, the same strategy applied to a chain of N 
identical quadrupoles in line (with transfer matrix T^), 
fails to detect any additional eigenfrequency. To under- 
stand this failure, we recall that eq. H18(l implicitly at- 
taches each of the two outgoing legs of the chain to the cor- 
responding incoming leg, thus enforcing periodic bound- 
ary conditions under which the system acquires an addi- 
tional translational symmetry along the chain. The low- 
est eigenmode of the system possesses the full transla- 



tional symmetry and shows the same current distribu- 
tion for each of the TV quadrupoles, implying an eigen- 
vector collinear to the current induced by the electrodes, 

Icxt (1, 1, ■ ■ • J 1)"^- Since all other eigenvectors are or- 
thogonal to the groundstate, they cannot be seen by sim- 
ply looking "if what goes in comes out" . (This fact is rig- 
orous for periodic boundary conditions; for closed bound- 
ary conditions, as assumed in Sec. the same mecha- 
nism gradually suppresses the response of the higher eigen- 
modes who have very little overlap with the groundstate 

— see Fig. 0) 

We will now present a new resonance condition which 
remedies these shortcomings; namely that, at resonance, 
all 4 eigenvalues /i of the transfer matrix of the entire 
chain should be unity. At this point, one may think of 
it as a necessary condition which prevents the norm of 
the state vector J from diverging or going to zero after 
many transfers through the same quadrupole chain. The 
proof that the new resonance condition is equivalent to the 
usual definition — i.e. a divergent impedance between the 
electrodes (and thus zero net current through the sample) 

— is deferred to the end of this section. 

In order to simplify the following calculations, we in- 
troduce the hcrmitian and unitary matrix 



(19) 



which allows for a transformation of the state vector to 
a symmetric basis, 3k = PJfc = {^k ^-^t ^^k) '^ith 
components U^ = -^{Uk ± Uk') and = -^{Ik ± Ik')- 
The associated transfer matrix of a single quadrupole in 
this basis is S = PTP, and thus of the general form S = 

A B 




( c D ) "^ith the 2x2 submatrices 



A 
C 



1 a 
1 

ci 

C2 



B 
D 



hi 62 


di d2 



(20a) 
(20b) 



Due to the particular structure of the submatrices — im- 
plying A" = ( J Y), AB = B, CA = C, and CB = — 
the transfer matrix of a chain of N identical quadrupoles 
can be calculated explicitly; as may be proved by induc- 
tion, it reads 



J TV 



.N 



N-1 

B E GfeC 
GatC 




(21) 



with the geometric series 



k-l 



Gfc ^ DP = (1 - D)-i (1 - D'^) (using Go = 0). 



p=0 



(22) 

(Obviously, the second equality from the left holds only 
as long as (1 — D) is invertible.) 
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has the same sparsity pattern, i.e. the same distri- 
bution of zeros, as S. In particular, its upper left 2x2 sub- 
matrix reads (qI), thus preserving the right eigenvector 

(1,0,0,0)"^, standing for the invariance to a global volt- 
age shift, and the left eigenvector (0,1,0,0), responsible 
for current conservation, both with eigenvalue fi ~ I. 

The remaining two eigenvalues are fi and 1/fi (since 
detS^ — det S = 1), and have their origin in the lower 
right submatrix, D^. Using = DtrD - IdetD (valid 
for any 2x2 matrix) , together with det D = 1 , we have 



parallel, Ta{^, 1), with ctl/ctc = -l/iLCco^) = 
another vertical capacitor, and finally two capacitors in 
parallel, Ta(1, 1): 

S = P Ta(1, 1) Tb(2) TA(f^, 1) Tb(2) P 





Vo 



, icj2 _2 + 2cj^ 
1 

19 - 10(^2 



5,~,2 



-16 




(26) 



N-1 ■nN-2 



(23a) 



A straightforward calculation shows that, in this case, the 
resonance condition H25(l reduces to 



with ^ = ^trD. The recurrence relation H23a|l may be 
thought of as a matrix version of the one defining the 
Chebyshev polynomials, eq. (|3J|. It may be easily verified 
that 

= Un-i{0^ - UN-2m , (23b) 

(with UniO the Chebyshev polynomials of the second 
kind, eq. ©) fulfils eq. (|23a(l and meets the initial condi- 
tions D° — 1 and = D. Quite similarly, for the trace 
of D^, (|23a|l gives 



/II 1 2TTm 

Urn — i\ cos 

6 6 N 



, / 5 1 . 9 rmr 
±W- + -sm^-, (27) 



trD^ = 2^trD'''-^ -trD 



(24a) 



which, along with the initial conditions trD" = trl = 2 
and trD^ = 2^, yields 



with m = 0. . .iV — 1. The only difference between the 
eigenfrequencies (|27|) . calculated for periodic boundary 
conditions in both, x- and y-direction, and the former 
result, eq. H12|l . obtained for a closed chain with peri- 
odic boundary conditions only in j/-direction, resides in 
the argument of the sine function whose symmetry about 
^ organises the eigenfrequencies (|27|) in degenerate pairs 
(except for uq and, for even N, iijv/2)- 

The second example to be analysed with the present 
method is the array of staircases shown in Fig. ^ and 
obtained from the unit cell in Fig. I^d) . For its transfer 
matrix. 



trD^ = 2T^(e = cos0) = 2cos[7V0] , (24b) S = P TA(f^, 1) Tb(1 + f^) Ta(1, f^) Tb(1 + f^) P, 

(28) 



where Tat are the Chebyshev polynomials of the first kind. 

Since the trace of a matrix is invariant under basis 
transformation, the resonance condition, /i = 1, may be 
reformulated as trD^ = ^+1/^ = 2, insertion of which 
in eq. (|24b|l requires 9^ = with m = . . . iV - 1, and 
thus finally: 



-trD 
2 



^ — cos 6j, 



27rm 
N 



the resonance condition H25|l is a fourth order equation in 

cZ;* - SCb^ + (18 - 4/3,2„) [b^ ~ 8^^ + 1 = 
[Cj^ - (4 + 2pm) + 1] [^^ - (4 - 2/3„) + 1] = , 
with f3m = ± cos ^ ; its solutions are 



(25) 



= A„ + 2 ± A/(/3m + 2)^-1 



(29) 



The resonance condition (|25|) is the main result of this 
section; it states that, in order to compute all resonances 
of a chain of N quadrupoles under periodic boundary con- 
ditions, one simply has to calculate the transfer matrix S 
of a single quadrupole — which, of course, depends on the 
components constituting the quadrupole and their setup 
— and then resolve the algebraic equation H25() for the 
lower right 2x2 submatrix D. 



3.1 Applications 

To see this recipe at work, let us recalculate the resonances 
of an array obtained by replication of the unit cell of 
Fig-EJc)- Using the capacitors' conductance as a reference, 
Co — <^c — i'-^C', the quadrupole's transfer matrix T may 
be written (from right to left) as a product of the transfer 
matrices of a vertical capacitor, Tb(2) (with twice the ca- 
pacity C, due to the periodic boundary conditions along 
the electrodes) , followed by an inductor and a capacitor in 



(where, again, m e [0;iV/2] suffices, since the rest of the 
resonance frequencies is obtained by symmetry). 

Fig. [3 displays the impedance of the staircase array of 
Fig. ^ for which = 3. A comparison of the resonance 
frequencies (|29|l (dotted lines) with an ENR calculation for 
the same array (continuous lines) shows that eq. H29() not 
only predicts the right number of resonances — this result 
is non-trivial since, without degeneracy, eq. (|29|l would 
produce 8 resonances ( 4 for m = and another 4 for m = 
1) instead of 7 — but also matches the position of most 
of the resonances. The deviations, observed in particular 
for UJ2+ and 1^3+, are due to the assumption of periodic 
boundary conditions along the chain, an approximation 
used in the derivation of eq. 



3.2 Link to the usual resonance condition 

Up to now, the resonance condition 125|l had to be thought 
of as a necessary one, since = 1 for all four eigenvalues 
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Fig. 7. Impedance between the electrodes for the setup 
shown in Fig. 01 with an array generated from the unit cell 
in Fig. Od), replicated 3 times in x-direction and twice in y- 
direction, embedded in a 9 by 4 capacitor lattice. The real 
part of the ENR impedance, calculated for L = C = 1 and 
Rcoii = 10~*, is plotted with continuous lines. The dotted lines 
indicate the resonance frequencies obtained analytically from 
eq. (1291 for = 3 and m = ((Do = 1 (doubly degenerate) 

and Cji± = \/3±2^), and m = 1 (ii2± = \j5/2 ± ^21/2 
and £l)3± = j3/2±V5/2 



For the groundstate, m = 0, the hne of reasoning is 
slightly more subtle: the naive procedure of solving the 
matrix equation (|18() at resonance, Ji + Ii' = 0, amounts 
in the present language to solving 

S»(!)^(!). (32, 

The right eigenvector {U, 0, 0, 0)""" may be subtracted from 
the system of equations, and one is left with an eigenvalue 
problem for the last two components of J, ( y ) = ( ^ ), 
which — using eq. — reduces to finding the right 
eigenvector of the single quadrupole's D, associated with 
eigenvalue fi — 1: 

D(^) = (^) (33) 

(Note that, for the generally asymmetric matrix D, the 
right eigenspace contains only one eigenvector if — as in 
the present case of a resonance — the eigenvalues are de- 
generate.) 

In all cases — for the groundstate and for the ex- 
cited modes — the condition that shall only have 
eigenvalues fj, — 1 is equivalent to the usual definition of 
a divergent impedance between the electrodes, implying 
Ik + Ik' =0 at finite applied voltage. 



of was merely required for the norm of the eigenvec- 
tors J to remain finite after a great number of transfers 
through the same chain of N quadrupoles. The aim of this 
paragraph is to show the equivalence of the new resonance 
condition (|25|l with the more familiar one, namely a van- 
ishing total current through any cross-section of the chain, 
I^ = -^{Ik + Jfe/) = 0, despite finite apphed voltage. 

At resonance, may be evaluated explicitly: from 
eq. (E5|l . 0rn = the recurrence relation l|23b(l . and 

the Chebyshev polynomials (jB)) evaluated at resonance. 



f/Ar_i(cos 9m) 

Un-2{cOS 9m) 



- [N9„ 



N if m = 0, 
else 



sin [{N - 1)6'„ 



N -I if TO = 0, 
— 1 else 




^jv )NT> - [N -1)1 if m = 0, 



else (i.e. m = 1 . . . — 1). 

(30) 

For the excited eigenmodes, to mod iV 7^ 0, = 1 
implies G„i = (eq. 122|) '). and hence, from eq. (|21|l . 



1 



(31) 



Generally a is finite; it is thus obvious that any vector 
J is a right eigenvector of with eigenvalue /i = 1, if 
(and only if) its second component, describing the net 
current through the sample, vanishes. 



4 Clusters in an infinite lattice 

The systems studied in the preceding sections could be 
reduced to effective ID problems, because (i) the arrays 
covered the whole width of the electrodes and (ii) periodic 
boundary conditions were assumed in y-direction. In this 
section, by contrast, we will examine simple regular ID 
arrays in an infinite square lattice. The observed changes 
turn out to be considerable, and sometimes not only shift 
the frequencies, but qualitatively change the spectrum. 

The first example we want to inspect consists in a sin- 
gle line of alternating inductors and capacitors embedded 
in an otherwise pure capacitor lattice, as shown in the left- 
most graph of Fig. |H1 Alternatively, this line is generated 
by replicating the unit cell of Fig.[5fc) only in a::-direction. 

Its eigenfrequencies are obtained very accurately with 
the so-called spectral method, initially proposed by Stra- 
ley Bergman ^U] and Milton j^, and later adapted 
to 2D lattices by Clerc, Giraud, Luck and co-workers jl2l 
113) . employing a Green's function formalism developped 
by McCrea and Whipple in the context of random 
walks; (see also Spitzer ^S))- This approach connects the 
impedance resonances of a (not necessarily ordered) clus- 
ter to the eigenvalues of a non-symmetric square matrix, 
which — in general — have to be evaluated numerically. 
The latter version of the spectral method assumes the 
cluster to be located in the middle of an infinite lattice 
— just as in the examples we wish to study. To obtain 
similar results from an ENR calculation one would have 
to incorporate the cluster in a much larger capacitor net- 
work (compared to the cluster's size) — a situation which 
is computationally very demanding. 
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Fig. 8. Low-frequency renormalisation scheme for the transformation of an infinite line of alternating coils and capacitors in a 
2D capacitor lattice with bonds of capacity C (left graph) to an infinite ladder of capacitors with only one horizontal capacitor 
replaced by a coil (right graph). (In the second and third graph, a bond with two capacitors in series always represents two 
times 2C in series.) 
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Fig. 9. Eigenfrequencies uj in units of luo for a line of N alter- 
nating coils and capacitors embedded in an infinite capacitor 
lattice. The circles represent numerically calculated eigenfre- 
quencies, obtained with the spectral method. The dashed lines 
are the analytical large asymptotes, (Ijlf 
0.7395 and Qhf =i l/VVS - 



l/\/2V2- 1 ~ 
1.1688 discussed in the text. 



The eigenfrequencies of the line of alternating coils and 
capacitors, obtained via the spectral method, are plotted 
with dots in Fig. IHl as a function of the number of coils, 
N. In the limit of large N, the lowest eigenmode shows 
the full translational invariance of the system, and thus 
has the same current distribution in every vertical stripe. 
This implies that the left and right boundary of the stripe 
can be thought of as attached together. In order to per- 
form this operation, illustrated in Fig. |S1 one has to (i) 
chose a stripe which is symmetric about one of the coils, 
(ii) replace all horizontal capacitors of capacity C in all 
pure capacitor columns by two capacitors in series, with 
capacity 2C each, (iii) cut between the doubled capac- 
itors, and (iv) tie the corresponding vacant ends of the 
stripe together. One ends up with an infinite ladder with 
horizontal steps of 2C (two times 2C in series, the whole 
parallel to a capacitor C), except for one step which con- 
tains a coil instead of the single capacitor; the vertical legs 
of the ladder consist of capacitors C {cf. rightmost graph 



in Fig. inj. The capacities may then be summed up in a 
procedure analogous to eqs. ((T^ and ((TH) . yielding a to- 
tal capacity of C(2\/2— 1), and thus the low-frequency 

asymptote cDlf — l/\/2\/2 — 1 ~ 0.7395 (lower dashed 
line in Fig.EJ. 

The position of the high-frequency asymptote cjhf — 

1/VV3 - 1 ~ 1.1688 (upper dashed line in Fig.©, on 
the opposite side of the spectrum, can be calculated al- 
most analogously: in the limit of an infinite LC chain in 
its highest eigenmode, the currents through two neigh- 
bouring coils are at any moment antiparallel, but of the 
same magnitude. The capacitors separating the coils hor- 
izontally are thus located at current nodes and can be 
omitted. The corresponding ladder looks the same as the 
one shown in the right of Fig.|Hl except for the fact that all 
capacitors of 2C (on the bent lines) have to be removed. 

Comparison to the case with periodic boundaries in y 
direction, studied in Sec.[21(eq. H12|l in particular), shows 
that the shift of a resonance depends on its location in 
the spectrum: at the low-frequency edge, formerly located 
at UJ ~ \/5/3 ~ 1.291, now at (D ~ 0.7395, the renormal- 
isation is stronger than at the high-frequency boundary 
(formerly tD ~ a/2 ~ 1.414, now ui ~ 1.1688). The rea- 
son for this behaviour lies in the different renormalisation 
schemes which — as pointed out above — couple the in- 
ductors to an equivalent capacity which is less enhanced 
for high frequencies than for low ones. 

The second cluster we are going to study in this sec- 
tion is a straight line of N coils in series embedded in 
a pure capacitor environment. The resonance frequencies 
of this "lattice worm" , obtained via the spectral method, 
are plotted with dots in Fig. ^1 For large N and at high 
frequencies, the cluster can again be reduced to almost 
independent stripes. The corresponding renormalisation 
scheme is illustrated in Fig. 1111 and ends up in an infinite 
ladder with legs of capacity C/2 and steps of capacity C 
(except for the central step which is a coil of inductance 
L). Summing up each semi- infinite capacitor ladder yields 
Ci„f = iC(l + V2) {cf. eqs (O and lO), and thus the 

high-frequency asymptote (Dhf — vl + \/2 ~ 1.5538 (up- 
per dashed line in Fig. Il()|l. Comparison to the periodic 
case, eq. where cu \/3 for to = and N 00, 
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Fig. 11. Higli- frequency renormalisation scheme for an infinite line of coils in a 2D capacitor lattice (left graph), ending up in 
an infinite capacitor ladder with only one coil (right graph). 
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Fig. 10. Eigenfrequencies uj in units of ujq for a continuous 
line of A'' coils embedded in an infinite capacitor lattice. The 
circles represent numerically calculated eigenfrequencies, ob- 
tained with the spectral method. The dashed lines are the large 

A'" asymptotes, iIjlf — \/l.l6/N and iDhf — + ^2 ~ 1.5538 
(see text). 



shows that the renormaUsation shifts arc rather modest 
at high frequencies. 

At low frequencies, by contrast, the infinite environ- 
ment has drastic effects: as can be seen from Fig. EH the 
lowest resonance tends to zero frequency as the length of 
the worm increases. The corresponding density of states 
therefore shows no gap at low frequencies — as opposed 
to the density of states (|10|l of the periodic system, illus- 
trated in the last graph of Fig. El 

The reason for this behaviour can be qualitatively un- 
derstood in a simple picture representing the whole worm 
of N coils as a dipole of total inductance NL. This dipole 
is thought to be coupled to a single capacitor bearing the 
entire capacity of the lattice. In this model, the lowest 
resonance would thus be found at w ~ y/ 1 / {N LCut) ■ As 
pointed out in Appendix ^ the dipole model evaluates 
the entire lattice capacity to roughly Ciat ~ C, a value 
which allows to reproduce the exact resonance frequency 
(D = 1 for a single coil {N = 1) in a 2D capacitor lattice. 
For increasing A'', Cut is found slightly reduced since more 
and more capacitors are replaced by coils. In the limit of 



large A^, the numerically calculated groundstate fits very 
accurately to uj ~ y/lAG/N (lower dashed line in Fig. llO() . 

In the same manner, the dipole model predicts a reso- 
nance frequency of tOm oc y^m/N for the m.th mode, with 
a proportionality constant of the order of unity, since the 
current nodes split the "worm" into m almost independent 
segments of N/m coils each. Comparison to the numeri- 
cally evaluated eigenfrequencies shows that this scenario 
works well for m < 5, but turns out to be too simplistic 
for the higher excited modes. 

On more general grounds, and observing that second- 
order terms do not improve the reproduction of the data, 
one may try an ansatz of the form 



(34) 



A linear dependence of the parameters, am ~ —0.40 + 
1.57m and Pm = -0.56 + 0.98m, deduced by a fit for A^ = 
200 . . . 1000 and m = 1 . . . 30, reproduces the numerically 
calculated resonance frequencies very accurately for a wide 
range of A^ and m, as long as N/m is large (> 5). 

We also note that the fitted a, with a slope of 1.57, 
concords with an analytical result obtained by Clerc et al. 
|12j . stating that Cjm ^ y/^m/N for the bulk states of a 
long linear cluster, i.e. if both m and A" are large. 



5 Conclusions and outlook 

In this article, we have studied the dielectric resonance 
spectra of ordered passive arrays, typically — although 
not necessarily — constituted of inductive and capacitive 
elements. Similar arrays, based on split-ring resonators, 
have recently been used to assemble metamaterials for 
the microwave regime, exhibiting negative refraction and 
other exotic properties 

In the first part, we have calculated the resonance fre- 
quencies of several arrays by solving a system of differen- 
tial equations deduced directly from Kirchhoff's rules — 
a technique which allows for a straightforward interpreta- 
tion of the spectra and thus provides a good handle for 
the influence of each of the array's parameters. On the 
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other hand, each new circuit requires a separate individ- 
ual analysis, relying on the solution of increasingly com- 
plicated recurrence relations as the unit cell of the array 
gets richer in structure. 

An alternative approach, presented in Sec. 13 deduces 
the resonance frequencies from the array's transfer ma- 
trix, i.e. a matrix connecting the state vectors (a combi- 
nation of currents and voltages at each vertex) at both 
ends of the cluster to each other. Within this formalism, 
the array is shown to be in resonance if all eigenvalues of 
its transfer matrix are unity for a given frequency — a 
condition which is demonstrated to be equivalent to the 
more familiar definition of a divergent impedance for a 
cluster between two electrodes, namely vanishing net cur- 
rent through the sample at finite applied voltage. Even 
large arrays, with complex unit cells, can be easily anal- 
ysed with this algorithm since, in a handier reformulation, 
the new resonance condition, eq. (|25(l . does not require the 
computation of the transfer matrix of the whole array, but 
only of a single unit cell. The latter is most conveniently 
evaluated by multiplication of the transfer matrices of a 
few standard situations — two in our case, cf. eq. (|17all . 
where the 2D arrays reduce to two-legged ladders due to 
the assumption of periodic boundary conditions along the 
electrodes — a task which can be performed using sym- 
bolic computation software. In its final form, the resonance 
frequencies are typically given as the roots of a polynomial 
whose degree depends on the complexity of the unit cell. 

If, however, the cluster happens to be embedded in a 
much larger lattice, periodic boundary conditions are not 
a pertinent approximation, and the resonance frequencies 
are generically not available in closed form anymore. In 
this case, major changes in the spectrum have to be ex- 
pected, of which some can be understood in terms of renor- 
malisation, while others cannot. Among the latter, one 
may recall the example of a longer growing linear array: 
embedded in a periodic medium, the spectrum remains 
gapped for any length of the cluster; if, by contrast, the 
same cluster is isolated in a homogeneous network, the gap 
closes with increasing length and the resonance spectrum 
carries the signature of one or more rather independent 
two-dimensional dipoles. 

The tools discussed in this article may be useful for the 
design of circuits with a custom electric response. Such 
devices could be used for labelling — just as bar codes — 
which could be read by an automated system. 

For the future it would be undoubtedly desirable, es- 
pecially in the perspective of the already mentioned devel- 
opment of metamaterials with negative refraction, to take 
into account the magnetic part of the optical response, 
with the aim to create devices having both, tunable di- 
electric and magnetic resonances. 
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A Dipole approximation for linear clusters 

The scenario presented in the following is based on the 
idea that the current flowing through a linear cluster iso- 
lated in a capacitor lattice flows back through the lat- 
tice with a current distribution which, in the lowest eigen- 
mode, resembles the field lines of a dipole in 2D electro- 
statics (see Fig. lT^ . 




x=0 



X 



Fig. 12. Current distribution (dashed lines) of a linear cluster 
of length A'^ = 1 in the dipole approximation. (The lower half- 
plane has been omitted for clarity.) 



This idea is supported by the spectral method, in the 
framework of which the electric potential induced by 
a single horizontal coil, located at the origin, is known to 
be given by ^2] 



-W G 



(35) 



where A is a frequency-dependent proportionality constant, 
W the voltage measured between the ends of the coil, and 
Gr the 2D lattice Green's function. Far from the coil, for 
|r| ^ 1, the lattice Green's function may be replaced by 
its continuous counterpart, Gr — 
reduces to a 2D dipole potential 



--^ln|r|, and eq. 



^W^G. 
ox 



W 



2tt + 



W cos ip 
2^ Irl 



(36) 



In this approximation, the current distribution follows 
the field lines ^ 

r((^) = — sin(/3 (37) 

TT 

where I is the length of the field line (in units of the lattice 
spacing). In order to calculate the resonance frequency of a 
cluster embedded in an infinite capacitor lattice, one has 
to assign an equivalent capacity, Ciat, to the entire lat- 
tice. Of course, Ciat depends on the current distribution: 
in general, capacitors almost perpendicular to the current 
lines contribute only little to Ciat, while capacitors follow- 
ing the current lines contribute as if they were in series. 
At great distance from the cluster, the lattice behaves as 
if it were assembled by independent lines of capacitors. 
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hence all in parallel, contributing each as the inverse of its 
length: 

k ^ 

The sum in eq. (|38|l runs over all possible current 
paths, and — in order to avoid divergences — the count- 
ing has to be done very carefully. In our example of a 
single coil in an infinite lattice, this may be achieved by 
supposing that the current can only leave the cc-axis at 
spots where a vertical capacitor is present. Following the 
illustration of Fig. 1121 this amounts to retaining only cur- 
rent lines passing through {x = k -\- = 1) with k = 
0, 1,2,.. .. In coordinates of the field lines H37|) . these points 
are described by 



h = 



1 



(39a) 
(39b) 



The first equation can be used to select only the desired 
current lines, as shown in Fig. 1121 Pinning the field lines 
in this way has a twofold virtue, namely on the short side 
to set a lower boundary for the path length l^. in H38|l . and 
to avoid the logarithmic divergence on the long side, Ik 
oo. After substitution in H38|l . and taking into account 
contributions from the lower half-plane, we obtain 
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TT J l + P 




= C 



(40) 



(The sum in the first line may be thought of as a midpoint 
trapeze approximation for the integral in the second line) . 

According to the dipole approximation, an isolated 
coil in a 2D cap acit or network thus resonates a,t lo — 



l/v^Ciat — 1/VLC = loq — which is the exact result 
known from the spectral method. 
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